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Abstract 

We construct cosmological spacetimes with null Kasner-like singularities as purely 
gravitational solutions with no other background fields turned on. These can be recast 
as anisotropic plane- wave spacetimes by coordinate transformations. We analyse string 
quantization to find the spectrum of string modes in these backgrounds. The classi- 
cal string modes can be solved for exactly in these time-dependent backgrounds, which 
enables a detailed study of the near singularity string spectrum, (time-dependent) oscilla- 
tor masses and wavefunctions. We find that for low lying string modes (finite oscillation 
number), the classical near-singularity string mode functions are non-divergent for various 
families of singularities. Furthermore, for any infinitesimal regularization of the vicinity 
of the singularity, we find a tower of string modes of ultra-high oscillation number which 
propagate essentially freely in the background. The resulting picture suggests that string 
interactions are non-negligible near the singularity. 
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1 Introduction 

Understanding cosmological singularities in string theory is an important goal, and has been 
the subject of several investigations e.g. ^^M^UMU\BMMTilTMT5lMW^MTi\ 

Our work in this paper has been in part motivated by investigations [23l 121] involving 
generalizations of AdS/CFT where the bulk contains null or spacelike cosmological singularities, 
with a nontrivial dilaton e* that vanishes at the location of the cosmological singularity, the 
curvatures behaving as Rmn ~ 8m^9n^- The gauge theory duals are A/=4 Super Yang- 
Mills theories with a time-dependent gauge coupling g\ M = e , and [23J El] describe aspects 
of the dual descriptions of the bulk cosmological singularities. From the bulk point of view, 
supergravity breaks down and possible resolutions of the cosmological singularity stem from 
stringy effects. Indeed noting a' ~ 3 1 N from the usual AdS/CFT dictionary and extrapolating 
naively to these time- dependent cases with a nontrivial dilaton, we have a' ~ indicating 
vanishing effective tension for stringy excitations, when e* — > near the singularity. While 
this is perhaps wrong in detail, we expect that stringy effects are becoming important near the 
bulk singularity, corresponding to possible gauge coupling effects in the dual gauge theory. It 
is therefore interesting to understand worldsheet string effects in the vicinity of the singularity. 
Owing to the technical difficulties with string quantization in an AdS background with RR 
flux, we would like to look for simpler, purely gravitational backgrounds as toy models whose 
singularity structure shares some essential features with the backgrounds in the AdS/CFT 
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investigations. We first find such spacetime backgrounds as "near-singularity" solutions to 
Type II supergravity (preserving a fraction of lightcone supersymmetry) . In general, these 
are null Kasner-like solutions with null cosmological singularities (we also find approximate 
solutions that extrapolate from these near-singularity solutions to flat space asymptotically). 
These can be recast as anisotropic plane-wave-like spacetimes by a coordinate transformation, 
and we outline arguments in these coordinates suggesting the absence of higher derivative 
curvature corrections to these spacetimes (which are essentially plane- wave backgrounds). 

We then perform an analysis of string quantization to find the spectrum of string modes 
in these backgrounds. We find it convenient to use (Rosen-like) coordinates where the null 
cosmology interpretation is manifest. With these lightlike backgrounds, it is natural to use 
lightcone gauge. The classical string modes can be solved for exactly in these time-dependent 
backgrounds, which enables a detailed study of the near singularity string spectrum. For various 
families of singularities, the classical string oscillation amplitudes for low lying oscillation num- 
ber n are non-divergent near the singularity, with asymptotic time-dependence similar to the 
center-of-mass modes. From the Hamiltonian, we find time-dependent masses for these string 
oscillator modes. However, for any infinitesimal regularization of the vicinity of the singularity, 
say r < r e , we find string modes of ultra-high oscillation number n 3> which propagate 
essentially freely in the background. The near-singularity region thus appears to be filled with 
such highly stringy modes. There have been several investigations of string quantization in 
plane-wave backgrounds with singularities [HI El [HI EHl EH] and our string analysis has some 
overlap with [8] in particular. 

In sec. 2, we describe the spacetime backgrounds. Sec. 3 describes the string quantization. 
Sec. 4 contains a discussion and open questions. In Appendix A, we describe some properties 
of the spacetime backgrounds, while Appendix B outlines string quantization in coordinates 
corresponding to a different time parameter. 

2 The spacetime backgrounds 

We are interested in spacetime backgrounds that have a Big-Bang or Big-Crunch type of sin- 
gularity at some value of the lightlike time coordinate x + . We also want to restrict attention to 
purely gravitational solutions for simplicity, i.e. with unexcited dilaton and RR/NSNS fields. 
This means we want to solve the equations Rmn = 0. Null-time dependence reduces these 
equations to R ++ = 0. 

Let us begin by considering a spacetime background with two scale factors, of the form 

ds 2 = g^dx^dx" = e f{x+) (-2dx + dx~ + dx l dx l ) + e h(x+) dx m dx m , (1) 
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where i = 1,2, m = 3, . . . , D — 2. One may also think of the x m directions as compactified, 
representing say a T D ~ 4 . For the critical superstring with D = 10, we could alternatively 
replace this 6-dim transverse space by some Ricci flat space such as a Calabi-Yau 3-fold. The 
intuition here is that the time-dependence of the "internal" space induces time-dependence in 
the 4-dim spacetime as well, as in e.g. [30l l3T]. Another perspective is that the internal space 
scale factor is the analog of the dilaton in the AdS/CFT cosmological solutions of [22] 121], as 
we will elaborate on below. 

Simple classes of singularities in this system are obtained for spacetimes whose limiting form 
in the vicinity of x + = is null Kasner-likt 

ds 2 = (x + ) a (-2dx + dx~ + dx l dx l ) + {x + ) b dx m dx m , a > . (2) 

More generally, consider spacetimes of a general Kasner-like form 

ds 2 = {x + ) a (-2dx + dx~ + dxUx 1 ) + {x + ) bm dx m dx m , a > , (3) 

i.e. the individual internal dimensions x m evolve independently according to their Kasner ex- 
ponents b m appearing in the individual scale factors e hm ^ x+ ^ — > (x + ) bm as x + — > 0. 
The coordinate transformation x 1 = (x + )~ ai ^ 2 y I , where a/ = a, b m , gives 

{x+Y>W = (V) 2 - aidx+y * dyI + aW J^ +)2 • (4) 

x + 4{x + ) 2 

Then the metric ([3]) becomes of manifest plane-wave form 

a 2 aj(a + 1) \ x 2 



ds 2 = -2(x + ) a dx + dy- + 



T o (» 



(5) 



where we have redefined y~ = x~ + (^^+^+1 )■ F° r a i — a ,b m distinct, these are in general 
anisotropic plane-waves with singularities (after further redefining (x + ) a dx + = dX). In what 
follows, we will find it convenient to work in the (Rosen) coordinates ([2D, ([3]), where the null 
cosmology interpretation is manifest, but as we will see below, there are close parallels with 
various previous studies on plane-wave spacetimes with singularities, most notably [8] (see also 

ESI 122!). 

The spacetimes 02]) have nonvanishing Riemann curvature components (with e.g. f = -^+) 

R +i+i = \ {(ff - 2/") ^ = ^f^(* + r 2 , 

R +m+m = \ (2fh' m - (h'J 2 - 2/4) = K^ + 2-b) (x+)b _ 2 (6) 



1 A spacetime of the form ([2]) ([3]) with a < can be transformed by a change of coordinates to one with 

a > by redefining y + = -^r. This recasts = (x + )~l a l = (y + )' a ' and moves the singularity at x + — > oo in 

the spacetimes with a < to y + — 0. Thus it is sufficient to study spacetimes ([2]) ([3]) with a > 0. 
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For these spacetimes to be Ricci-flat solutions of the Einstein equations, the equation of 
motion = must hold, giving 

r ++ = \wf -r+l E(- 2 ^ - (^) 2 + 2/XJ = o 

m 

=> a 2 + 2a+^^(-^ + 26 m + 2a6 m ) = . (7) 

m 

This relates the various (null) Kasner-like exponents a, b m . The equation in terms of the general 
scale factors shows that the curvature for the 4D scale factor e' is sourced by those for the 
"internal" scale factors e hm : indeed the h m are the analogs of the dilaton scalar in the AdS/CFT 
cosmological context [23J, [23] where the corresponding equation was R+ + = |(<9 + $) 2 . That is, 
the kinetic terms (d + h m ) 2 (and related cross-terms) play the role of the dilaton in driving the 
singular behaviour of the 4D part of the spacetime. 

In what follows, we will specialize to the symmetric case here, i.e. all b m = b equal (and 
e hm = e h ). Then R ++ = simplifies to 

-{ff-f' + ^-^{-2h"-{h'f + 2f'h!) = a 2 + 2a+^^(-b 2 + 2b + 2ab) =0 . (8) 
2 2 2 

If b = a, this equation (assuming D > 2) simplifies to give the solutions b = a = 0, —2, in which 
case the Riemann curvature components are seen to identically vanish (the solution (—2, —2) 
can be shown to be flat space by the coordinate transformation to plane- wave form). Thus an 
interesting solution requires that the "internal" x m -space either grows or shrinks faster than 
the spatial part of the 4-dim cosmology. For any 6 ^ a, the equation of motion above is a 
quadratic in a that admits various solutions with 

2a = -2 - (D - A)b ± y/A + (D - A)(D - 2)b 2 . (9) 

Taking the positive radical, it can be seen that restricting a > for our solutions implies b > 2 
or b < 0. Furthermore a + 1 — b > if b < or \b\ < -jj^- 

Suppose we focus on finding solutions with a, b, being even integers, so that the metric 
allows unambiguous analytic continuation from x + < to x + > across the singularity. One 
may imagine that this is a coordinate-dependent choice of the time parameter x + and therefore 
not sacrosanct: however if we do take this choice, a, b being even integers seems natural. This is 
more restrictive: we must consider the above as Diophantine quadratic equations with solutions 
over integers, which are in general rarer. We then need to look for those b for which the radical 
above is integral. For the cases of obvious interest, i.e. the bosonic string (D = 26) and the 
superstring (D = 10), the radicals simplify to 2a/1 + 1326 2 and 2vT^TT26 2 ~ respectively. It is 
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then straightforward to check that 

(a, b) = (0, 2), (44, -2), (44, 92), (2068, -92) ... [D = 26] , 

(a, 6) = (0, 2), (12, -2), (12, 28), (180, -28), (180, 390), ... [D = 10] , (10) 

are solutions. Our analysis of these solutions in what follows will not depend on these detailed 
values though. 

These solutions can be thought of as the "near-singularity" limiting regions of more general 
spacetimes where the scale factors e^,e hm are not necessarily of power-law typd^l Since the 
various scale factors e-^,e hm , are related by the single equation of motion fl7|), a generic choice 
of & admits a solution to for the remaining scale factors e hm . For instance, with a single 
scale factor e hm = e h , taking & = tanh a (x + ), we can in principle solve for e h . In the limiting 
near- singularity region, we have already seen null-Kasner-like solutions with (jSj) relating the 
Kasner exponents. In the asymptotic region of large x + , it can be checked that 

e f = tanh a (x+) -> 1 - 2ae~ 2x+ , e h ~ const + -^— e ~^ + (n) 

D — 4 

is an approximate solution to (jHJ) (dropping the subleading nonlinear terms). 

We now make a few comments on the cosmological singularities in these spacetimes. No 
curvature invariants diverge due to the lightlike nature of this system, since no nontrivial 
contraction is nonzero. However there are diverging tidal forces for null geodesic congruences. 
Consider for instance a simple class of null geodesic congruences propagating solely along x + 
(at constant x~,x l ,x m ), with cross-section along the x l or x m directions. These are described 
by (r+ + = /' = ^ is the only nonzero Vt) 

^\ Ttl m(*4)^ +n j^)\ . (12) 



d\ 2 13 \d\ J \d\ J d\ 2 ++ \ dX 
This gives the affine parameter along these null geodesies 

/+ r (x + ) a+1 
dx + e fi - x ) = const. / dx + (x + ) a = const.- — . (13) 
J a + 1 

and the tangent vector 

t = dx=(^?\d± = ?-d + . (14) 



The relative acceleration of neighbouring geodesies in a null congruence can be calculated using 
the geodesic deviation equation giving 

a M = g MN R NC BDt c en B (15) 



2 These solutions also arise as certain Penrose limits starting with some cosmological spacetimes and adding 
a spectator dimension [8] (see also [32]). 
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where n = n B ds is the separation vector along a cross-section of the congruence. For our 
system, this gives 



n i — n iir> /t+\2„i _ «(fl+^K m _ mm p { c+\2^™ — ^fta + 2 b)n" 

a - g R +i+i {£ ) n - ^ %+ y a+2 , a - g R +m+rn (£ ) n - 4 ^ 2a+2 ■ ( 16 J 



The corresponding invariant acceleration norms are 



\„i\2 _ „i „i \.,m\2 „ „m„m /'1'7^ 

a = qua a ~ -. — , .„ — , a = s mm a a ~ -. — — — r— — — , li 

1 1 1)11 ( x +\3a+4 ' 1 1 nmm , +\-fe+4 a +4 ' V / 



So we see diverging tidal forces as x + — > for spacetimes satisfying the conditions (restricting 
to a > 0) 

£><4a + 4, a>0, (18) 

indicating a singularity^. Since tidal forces diverge (for a, b, satisfying both these conditions) 
along both the x l and the x m directions, the locus of the singularity is the 8-dim space spanned 
by the x\ x m . From the point of view of a Penrose-like diagram, we see that the singularity locus 
extends all the way to x~ — > oo. We will see reflections of this later in the string worldsheet 
analysis. 

In Appendix A, we show that these spacetime backgrounds preserve 16 real (lightcone) 
supercharges. This is not a feature we use however, and our worldsheet analysis below does 
not appear to depend crucially on spacetime supersymmetry of these backgrounds. 

We also mention that these spacetimes appear to not admit a' corrections due to higher 
order curvature terms, as is often the case with lightlike backgrounds. This is perhaps not 
surprising in light of the coordinate transformation that casts these null cosmologies in the 
form of anisotropic plane waves, which are known to be devoid of higher derivative corrections. 

In general, these spacetimes are slightly different from those studied by e.g.^ El IH El E] 
which were time-orbifold-like spacetimes. Although there are conceptual similarities, the de- 
tailed structure of the spacetimes are different and in particular, there is no issue of backreaction 
due to several orbifold "images" [5JE]- 

In what follows, we analyse the string spectrum in the vicinity of the cosmological singu- 
larities of these spacetimes. 



3 A string worldsheet analysis 

We will now describe a worldsheet analysis of string propagation in these backgrounds. Consider 
the worldsheet action for the closed string propagating in such backgrounds 

S = j drda V^hh ab d a X»d b X u g^(X) (19) 

3 This is true except when the coefficients of all a 1 vanish: this happens for the spacetimes (a, b) — (0, 0), (0, 2). 
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The worldsheet metric h a b has the signature (-1,1). It is convenient in the worldsheet analysis 
to use lightcone gauge x + = r, in keeping with the null structure of the spacetimes in question 
here. Unlike flat space however, it is not possible in general to use both lightcone gauge x + = t 
and conformal gauge h a b oc r) ab since that is one gauge condition too many, as we will see below. 
Let us therefore begin by setting h Ta = 0, to simplify the worldsheet action, as in [33] (see also 
' ^ Then the worldsheet Lagrangian becomes 



£ = -^7 J da (-EgudrX^rX- 1 + ^g IJ d a X I d a X J -2Eg + _d T X-^j 



(20) 



where we have defined E(t, a) = y~j^-- Since X~ is not dynamical, we can eliminate this and 
reduce the system to the physical transverse degrees of freedom. Now if E = 1 is allowed, then 
we have h TT = —h acr , which is equivalent to conformal gauge being compatible with lightcone 
gauge. However, since the momentum conjugate to X~ is p_ = ^ g+ 7 which is a r- independent 
constant, we have E = (setting p_ = — ^7 by a r-independent reparametrization invari- 

ance). Thus we see that conformal gauge is disallowed^ since g + _ ^ — 1. The action for our 
background simplifies to 

S = -^— t J d 2 o ((drX 1 ) 2 - e^^id^X') 2 + e h{T) - f{T) (d T X m j 2 - e h{T)+f{r) (d a X m ) 2 ) (21) 

This action contains only the physical transverse oscillation modes X 1 = X 1 , X m , of the string. 
In effect, all the gauge freedom and corresponding constraints have been used up, with X~ = 

X Q +P-T. 

The corresponding Hamiltonian, — p + , re-expressing the momenta Il / in terms of d T X ! , is 

H = — *— / da \(d T X 1 ) 2 + e 2/(T) (d CT JT) 2 + e h{T) - f{T \d T X m ) 2 + e h{r)+f( - T) (d a X m ) 2 } (22) 
Ana.' J 

In general, one might imagine that a time-dependent background pumps in energy and excites 
string modes, and the classical Hamiltonian above does reflect this. For spacetimes satisfying 
e 2 ? — > near the singularity r — > 0, the potential energy of the X 1 modes due to the e 2 ^^ 
factor becomes vanishingly small near x + = (for a > 0). This could be taken to mean that it 
costs vanishingly little energy to create long strings as we approach x + = r = 0, the effective 
tension of string modes becoming vanishingly small near the singularity. However this appears 
to be misleading: what is relevant is e.g. the ratio 6 t^j^t ■ This has a more detailed form 
involving nontrivial r-dependence stemming from both gjj and from the asymptotic behaviour 



Ref. [35] studies some aspects of string quantization in Brinkman coordinates. 

5 Appendix B contains a discussion with the affine parameter A being the time parameter: in this case, 
= — 1, and conformal gauge is compatible with lightcone gauge. 
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of string modes X 1 , which we can solve for exactly in this background. Furthermore since this 
Hamiltonian corresponds to ^-translations and g + - ^ — 1, the string oscillator masses, which 
are coordinate invariant, are m 2 ~ g + P+P~, whose r-dependence is different from that of the 
Hamiltonian. In the case of affine parameter quantization (Appendix B), the time-dependence 
of the Hamiltonian translates directly to that of the oscillator masses. 

Heuristically one might imagine that the string gets highly excited and breaks up into bits 
propagating independently near the singularity: in a sense, this is akin to a worldsheet analog 
of the observations of BKL [36] on ultralocality near a cosmological singularity. It would be 
interesting to understand this better. We will find some parallels with this in our analysis later, 
which will reveal distinctly stringy behaviour. 

In the next section, we will study quantum string propagation in this background in detail. 
We will focus on the symmetric case, i.e. all b m = b equal giving two exponents a, b, but it is 
straightforward to generalize our analysis to the general case. 

3.1 String modes and quantization 

We are interested in studying the behaviour of string modes as we approach the singularity 
from the past, i.e. r < 0. For notational convenience, we will simply use r to denote |t| = —r 
in the expressions below. The equations of motion from the worldsheet action above are 

g2 x i _ e 2f(r) d 2 xl = q ? 

d 2 T X m + {d T h - d T f)d T X m - e 2fiT) d 2 a X m = , (23) 
which simplify in the near-singularity region of spacetime to 

d 2 x l - r^dlx 1 = o , 

d 2 X m + Q^m _ T 2aQ2 X m = q (34) 

T 

Decomposing the X 1 as /^(T)e m<T , we can show that the time-dependent mode solutions of 
these equations are given in terms of arbitrary linear combinations of two Bessel functions^ 

nT a+i\ — / nr a+i 



flW) = c',v/nr J 1 + c\o\/rvr Y 1 , 



(i 



Setting /* — * \frP n , /™ — > r v 'f™ ', transforms the equations of motion to the standard Bessel forms 



S 



These expressions are valid for v > 0, while similar Bessel functional forms with index J^y hold 
for v < 0. The Bessel index in /* is thus always less than ~ since a > 0, while for b < 0, 
the Bessel index in f™ is always greater than |. The complex coefficients c^ 1; c^ 2 can be taken 
to indicate the choice of a vacuum by defining positive/negative frequency modes. For now, 
we keep them as two independent unfixed constants: we will comment on specific choices at 
appropriate points in what follows. 

Note the similarity between these string worldsheet mode solutions and the well-known 
Hankel function description of spacetime scalar modes propagating in 4-dim de Sitter back- 
grounds. Spacetime scalar modes in the present null Kasner-like backgrounds are somewhat 
different however from thes^]. 

We can also examine the behaviour of the zero modes or center-of-mass modes. For n — 0, 
the equations of motion (j24p for Xq(t) can be solved to give 



X*(t) = + V2^a' Pl0 r , X m (r) = + ^a' Vma r 2v , (26) 

V 27T V Z7V 



where p/o are the center-of-mass momenta defined later (1291) . These show that for singularities 
with 2v > 0, the center of mass of the string is not driven to infinity by the singularity. We 
will find parallels of this with the asymptotics of low lying string oscillation modes. This is 
to be contrasted with the divergences for spacetimes with 2v < 0. Note that the zero mode 
behaviour is essentially point-particle-like. Thus the centers-of-mass of say a collection of 
infalling strings would appear to exhibit diverging tidal forces through geodesic deviation. 
However the crucial point is that the oscillations of the string are now non-negligible (even 
if finite). Thus neighbouring strings would appear to have large spatial overlap and string 
interactions become important near the singularity. 

The mode expansion for the spacetime coordinates of the string is 

oo 

71=1 

The constant will be fixed by demanding canonical commutation relations for the creation- 



7 Consider a massive scalar <\> in the background ([2]), with action S = J d D ' x^J — g(— g^ u d^(f>d u (j) — m 2 <f> 2 ) , and 
equation of motion —^=d tl (y^gg ,J '' J d v cj)) — m 2 (/) = 0. Taking modes 4> — f(x + )e tk - x + lk i x '+ lk ^ xTn j this simplifies 

t0 J-£+=2t ( k f + k m( x+ ) a ~ b + m 2 {x+) a + 2a+ ^T 4)b ) , which can be solved to give 

= exp [— + k^-^ + m 2 ^- + log, 

Thus generically these modes have a phase that oscillates "wildly" near the singularity x + — > 0. 
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annihilation operators. The momentum conjugates II 7 = d(/x 7 ) are 

= sb**. n-(r, = (28) 
We define the center-of-mass momenta p/o as 

Pio = / -7=^ = -^-XHr) , p m0 = -=IT = -^-X™(t) . (29) 

JO V27T \J2llOt' Jo \J 1-R yZTTOt' 

Then we see that imposing the nonzero commutation relations 

[X^pjo] = id j , [«n> a -m] = n b IJ Snm , [2n> 2 -m] = n ^ J &nm , (30) 

implies the equal time commutation relations, e.g. 

[X t (t, a), n J (r, a')) = ^-5 IJ 1 1 + ^( e m ^ CT ') + e -«(^)) J = ^ /J 5(a - a') , (31) 



n=l 



using the Fourier series representation for the Dirac (5-function, with the constant k T n being (this 
agrees with the conventions of [37] for flat space, except for a reversal of left /right movers) 



K n ~ \\ n\ I \l il^' ° n0 ~ C nl C n2 C nl C n2 ; l^J 

n y / l c n.oll a + -U 

where c^ is the Wronskian. We have used above the expressions for the derivatives of the mode 
functions f£ and some recursion relations for the Bessel functional to calculate the Wronskian 
of flfl 

Let us now discuss level matching. The operator that generates a-translations is worldsheet 
momentum P given by the stress tensor 



T ab ~ --1= ^ ~ - (^gjjdaX'd.X' - h ab h^gud^dixA 



(33) 



i We have used the following, the Bessel function argument being ( "L_ x ), 



/ \ 2i/ . . dJ v (z) 

J v -i{z) + Ju+i(z) = — J v (z) , J v -i{z) - Jv+i(z) = 2 ^ , 

Y v -i(z) + Y v+1 (z) = —Y v (z) , Y v _ x (z) - Y v+l {z) = , J v (z)Y v ^(z) - J v - X {z)Y v {z) = — 

z az irz 



10 



Then the a-translation gauge invariance is fixed by demanding that the momentum operator 
vanishes on the physical states, i.e. P = J daT Tr7 = 0. From our action above and our lightcone 
gauge condition h ra = 0, we have 

P = J da^d^d.X 1 + r b d T X m d a X m ) . (34) 

Using the mode expansion (1271) . this can be evaluated as 

P ~ r a J> (VlX - aLX) + (a™ a™ - a-„C)) , (35) 



n 



where we have used the Bessel recursion relations and the expressions for f£ (suppressing some 
overall unimportant numerical factors). This recovers the level matching conditions N = N. 
Now we calculate the string Hamiltonian. Using the mode expansion ( 1271) . we first evaluate 

n 

-{fnYW-n^-S) , 

i- ^ da{d <T X I f = ^n 2 |^| 2 (|/i| 2 ({^,aiJ + K,aU)-(/i) 2 K,^} 

n 

- (/«f{«- ni «-n}) • (36) 

The Hamiltonian (1221) then simplifies to 
H = ^((^) 2 + r^(X m ) 2 ) 

+ ^7 E l^l 2 ( (K, «-J + VJ) (\f n \ 2 + n\^\f n 

n \ 

- R, <} (tin) 2 + n 2 T 2a (f n f) - {ai n , ~at n } {{r:f + nV(tf ) 2 ) 

n \ 

" {CO (r 6 - a (/D 2 + ^^ +a (/D 2 ) 

- {a? n , a-} (r^UTf + nV +a (/r) 2 )) • (37) 
In the next section, we will examine free string behaviour in the vicinity of the singularity. 



(i 1 2 
'nl 
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3.2 Strings in the near singularity region 

Let us now understand the behaviour of the string mode functions near the singularity. It 
turns out that the near singularity limit r — > must be taken with care. We define a cutoff 
r = t € ~ as a short time regulator in the vicinity of the singularity r = 0. Then define 
n e = -^qrr as a cutoff on the worldsheet oscillation number. We then see sharp differences 
between the behaviour near r = r e of string modes with "low lying" oscillation numbers n < n e 
{i.e. nr^ +1 <C 1), and highly oscillating string modes with n 3> n e (i.e. nr^ +1 ^> 1). 

Noting the asymptotics J± u (x) ~ x ±v for x ~ 0, and Y u = cot(nv) J u — cosec(7r^) J_ u , we 
see that, near r = 0, the /^(r) approach, 

fn -> Ko + Kr r i fn -> Ko + Kt t2v ( r -> °) , ( 38 ) 

for modes with low lying oscillation numbers n < n e . The constant coefficients are (from the 
asymptotic Bessel expressions) 

n \^ <i + < 2 cot^ / n cosec^ 



v "\2a + 2j r(i±|) ' ™ ^ , \2a + 2) r(|±±) ' 

,m I „m „„f I^7T / \ „„„„„ VTT 



n V +1 C + Ccot^ / n \ q +1 cosec a+1 

2a + 2j r(2±^i) ' n0 n2V V2a + 2; r(*g^) 



Thus we see that the asymptotic r-dependence of such finite n string oscillation modes near 
t — > is essentially the same as for the the center-of-mass modes of the string fl26|) . Thus 
the (classical) string mode amplitudes are non-divergent near the singularity for cosmological 
solutions with 2v = a + 1 — b > 0. The string oscillation amplitude in such a curved spacetime 
is perhaps better defined as gmmifn) 2 '- this gives the asymptotics to be non-divergent for 

vl \I 



2a + 2 > b. In what follows, we will find useful the Wronskian combinations for the A^ , A^ 



cosec 7 * 



A* = \ l \ l * — \ l \ l * — nr l '2a+2 

V2a+2^ V2a+2^ 
cosec JdZL 

Am _ \ra \m* _ \m \m* _ „_m a+1 J _ J I* _ J* J ( A(\\ 

1 V a+1 J V a+1 / 

On the other hand, consider now modes with n 3> n e . Then we can see from the Bessel 
mode functions (|25|) (or directly from the equations of motion $M§) that these are oscillatory 
near the singularity: the argument n ^ cannot be taken to be small and the asymptotics (|38|) 
above are not valid. For instance, choosing linear combinations c^c^ = l,=t*, gives modes 
that are the analogs of ingoing or outgoing plane waves, i.e. the f£ are Hankel functions dressed 
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with powers of r, with asymptoticgj for r — > 

P L_ pdbtf»T»+V(a+l) ftn L p ±mr«+V(a+l) ^ ( 41 \ 

Note that for any regulator r e , however small, in the vicinity of the singularity, there exist 
modes of sufficiently high oscillation n such that the corresponding modes are of this form 
f |4TT) . Since the string oscillation number n can be arbitrarily large, such modes exist uniformly 
for all singularities, with 2v ^ 0, and are in a sense transplanckian: they are reminiscent of high 
frequency scalar modes propagating in an inflationary background. This behaviour, somewhat 
different from the finite n mode behaviour, is distinctly stringy. 

We first analyse the case 2u = a + 1 — b > 0. Using the asymptotic forms of the mode 
functions near the singularity r — > 0, for finite n < n e modes, 

fn Ko ! fn ~^ Kit j fn ~~ * Km j fn ~^ Kir ( 2z/ ) ^ 1 > (42) 

the Hamiltonian ( 1371) simplifies to 



+ ^ E \ k -\ 2 ( ((K. «- J + R> S!_J)|A 

n \ 



i |2 



KX}(AU 2 -{a^,5!_J(A») 2 ) 



+ n 2 r 2 « ((R,a l _ n } + R, alj) |A* | 2 - R,<}(A* ) 2 - R_ n , a*_J(A» ) 2 ) 
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+ E ( r a - 6 (2^) 2 ((«, a™ } + R, 3™ }) | A- | 2 - R, C)(A™ ) 2 - {a- , a^}(AT ) 2 ) 
+ n 2 r b+a (({Coy + {<,^J)ra 2 - WRKA™) 2 - R n ,^J(A™) 2 ) j . (43) 

Note that there are "interaction terms" of the form a^a^ and a^a 1 ^ besides the diagonal number- 
operator terms. The interaction terms have the same r-dependent coefficients as the diagonal 
terms so that they are not unimportant and cannot be ignored. 

The corresponding calculation for flat space (a, 6 = 0) involves sine and cosine modes (the 
analogs of the Bessel-J, Y), the Hamiltonian having no time-dependence. Analysing this near 
r — > 0, we see that the coefficients of the a^a^- and a^ n a^ n -terms are of the form (c 2 + c 2 ,) and 
(c* 2 + c^j 2 ), while that of the diagonal terms is (|ci| 2 + | C2 1 2 ) , where ci,C2 are the coefficients 
of the sine, cosine: then we see that choosing the usual positive frequency modes with c±,C2 



'This is also the asymptotic behaviour near r — * 00 of the modes (|25p for any n. 
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being 1, — i, results in just the diagonal term in the Hamiltonian. In the present case, due to 
the extra r-dependences in the Hamiltonian, the resulting expressions do not simplify and the 
"interaction" terms remain. A similar calculation with different choices of the basis modes (e.g. 
Hankel functions) yields equivalent results. 

This Hamiltonian ( )43l) . corresponding to the choice of time coordinate^!, can now 

be recast as 

H = W(fa ) 2 + r a - b (p m0 ) 2 )+ E 2(Q +i V (^(^+ n'r^Xo) 

+ JL ({2vfr a ~%^bZ + nV+Xcto) ) , (44) 

i c noi v 7 y 

where we have defined new oscillator modes (and their Hermitian conjugates) 

b L = A nO«n - A n0 «-n . b L = A nr«n ~ X nr O-n > 1 — 1,171 . (45) 

The string oscillator masses are Lorentz invariant expressions 

m 2 = -2g + -p +P ^ - g n {p I0 f . (46) 

From the above expressions, and recalling that p_ = — 2^7, ~P+ — H, we see that the center-of- 
mass terms cancel resulting in the time-dependent masses for these low-lying n < n e oscillation 
string modes 

m \r) = —L— T (1^- + T ^ + <>1^- + ^), [ 2 „>0], 
2a'(a+\) J- |cj T*n*|c&| |<S|/ ' - 1 ' 

(47) 

defining 

t-Ct, Ko = b%b i n0 , NZ = b^bZ, ACo = OCo- (48) 

These expressions should be understood as valid in the vicinity of the singularity, but only upto 
the regulator r < r e . 

The original left- and right-moving oscillator operators can be re-expressed in terms of b ! n as 
1 h 1 - X 1 * h 1 \ n 1 - 1 

I \ A nr u n0 A n0 u nr) > "n ~ .J, 

l n,0r ly n,0r 

and the level matching condition ( 1351) is recast as 

= £ niaLX - aLX) = £ -rf- ( 6 nXo - KA) ■ (50) 



< = -j^- (At bi - X J n l bl) , ~ai = -±- (At ^ - A^ 6ft) , I — i,m , (49) 



A 7 



10 



The affine parameter quantization, Appendix B, yields similar results as we describe here. 
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The commutation relations satisfied by the b J n0 , b J nT are 

\PmOi b n o\ = = [b mT , b^l] = [b m Q, b J nT ] , [&m0; ^nr] = n ^nflT$ ^ ' $mn = ~ [& nT , bJ Q ] , 

WL, Ko\ = -nAioX^ IJ S m n , [N£ T , bH] = -nAi^Xl , 

[N^Ni] = n^S^A'XXr + &Mo) . (51) 



using the left- and right-moving a, a-oscillator commutators (1301) . and the Wronskian combina- 
tions 0r s from f|40l) . Since the 6^ , b^ T , operators commute with their conjugates, the N^ , N* T 
operators do not have a number-operator-like interpretation on states annihilated by 6^ , b^ T . 
From the expression for the time-dependent masses, it is tempting to speculate that states that 
have e.g. vanishing (N^ T ) but nonzero (Nln) will become massless near the singularity r — > 0. 
However since N^ T , N^ , do not commutd_l|, these are generically not simultaneous eigenstates 
of N^o and N^ T , or energy eigenstates. If such a possibility can be validated for these 6 / -states, 
then the b l n0 , &™ T -oscillator states are light near the singularity while the b l nT , fr^-oscillator states 
are massive near r = x + —>■ 0, for the singularities with 1v > 0, b < (while for b > singu- 
larities, the ^-oscillator states are light and the 6^ T -states are massive). All these are light 
relative to the typical curvature scale however, as we will outline later. Some description of the 
^/-states is given in the next subsection: it would be interesting to develop ths further. 

Let us now consider the case 2v = a+1 — b < 0. Then the modes f™ behave near r — > as 
/™ -> A™ r 2v , while /™ -> A™ (2i/) r 2l/ . Thus A™ does not appear in the Hamiltonian (|37j) 
evaluated near r — * 0, which thus shows all a m ,a m -terms having identical asymptotics with 
time-dependence as r — > e.g. 

r b-a|ym|2 + n 2 T b+a|ym|2 _ ({2vf + n 2 T 2a+2 ) -> (52) 

in the coefficients. It is therefore not particularly insightful to recast a m , a m in terms of the b m - 
operators. The invariant oscillator masses thus grow as i and \ for the b l nT - and a m -oscillator 
states. The 6^ -states are light as before. 

Now let us consider the high oscillation modes with n ^> n e = -^+t: these have a uniform 
behaviour for both 2v ^ 0. Then using the asymptotics (HT1) for such modes (with c l nX = 
1, c^ 2 = —i, which are positive frequency), we see that 

/n - (-^ - -) — ^ , C ~ (-mr« - - ^ . (53) 



11 In terms of the original a, a-operators, this expression is 

WnOi N nT ] — n5 IJ [(XnQ^ir + ^n-r^no) ( a -n a n + ®n®-n) ~ ^nXnXn^n ~ 2A„q A^* a_ n a_, 
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This is very similar to the asymptotics of the modes (1251) at early times \t\ — > oo: however 
we are considering a different limit here, with large n, small r, and nr a+1 3> 1, so it is worth 
elaborating a little. In this limit, we calculate the expressions in (13T1) and express them as 

J_ (tpY + n 2 r 2a ( f) 2 ) ~ r a ( - + - ^ e -MnT*+V(«+i) 

J_ f T b-a ( im\2 , 2 6+a/ ^ af ^ , tb ^ „-2mr«+ 1 /(a+l) 

n2 ^ Un J + ™ T {Jn ) ) T ^^^y + ^+1) J 

_L + n 2 T 2a \P n \ 2 ) ~ 2r a , 1 (r b -1/ n m | 2 + nV+«|C| 2 ) ~ 2r« . (54) 

The expressions in the first two lines are vanishingly small relative to the ones in the third, and 
the Hamiltonian (1371) simplifies to 

Hn»n e ~ T a TTT ( a -n fl n + °Lrfln + ™) , (55) 

Li ~r -L 

J; n>n e 

as for free string propagation. The overall factor r a arises as before from the fact that we 
are using x + as time coordinate (with g + _ = — (x + ) a ). The oscillator masses for these highly 
stringy modes, using (jl6"j) . become 

-9+-H-,= E -4t(^ + ^» + ")» (56) 
a' ' a + 1 

n>n e 



m 2 (r) ~ 



as for free strings in flat space. The zero point energy has an ultraviolet completion as in that 
case. Thus these highly stringy modes exhibit essentially free propagation in these backgrounds. 
Comparing the mass v/^ of a typical single excitation state with the typical curvature scale 
set by the tidal forces \a % \, \ a m \, in this region, we have —rru2 ~ n-r3 " +4 n ^ m- 4a + 4 ^ 
Thus states satisfying <n< t3 ] +4 and <Cn< r 4 a +4-b are light relative to the local 
curvature scale. Similar comparisons for the n <n e states with the r-dependences T ±a and r ±b 
relative to the typical curvature scale hold if b < 2a + 2. 

For any finite, if infinitesimal, value of the near-singularity cutoff r e , such highly stringy 
modes exist, for oscillation number n 3> n t = -^+r , although naively removing the cutoff 
would suggest the absence of any such modes. 



3.3 Near singularity string states and wavefunctions 

We describe here some aspects of string states near the singularity using our discussion in the 
previous section, beginning with the low lying oscillation mode {/-states. 

Noting that the b 1 operators are complex linear combinations of the a 1 ,^ 1 ^, and recalling 
Bogolubov transformations, a {/-vacuum \<p) (annihilated by the b 1 ) would appear to be a 
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multi-particle state in terms of the original a 1 , a 1 operators, and vice- versa. Indeed we have 

<0|^Ago> = 5>|A^| 2 , P = 0,r, (57) 

n n 

where o£|0) = = c£|0). Similarly, defining |& ) = ^ |0), [&£> = ^ |0), |6 r ) = £ T |0>, |fet) = 
&„ r |0), it is straightforward to show that the lowest excited states have 

WLPmM) =0 ' ( b p\ b i q b m q \ b r) = U Kr Kp (^2 m \Knq\ 2 + n \Kiq\ 2 ) , P, 9, T = 0, T , (58) 

m 

using the expressions for the b 1 in terms of the a 1 , a 1 , and their commutation relations. 

Now it can be shown that [N^ , (b^ T ) 1 ] = /(nA^* 0r ) S IJ 5 nm (b 1 ^) 1 ' 1 b^ , using the ^-oscillator 
algebra fl5Tl) . Assuming the existence of a 6o _vacuum 7 defining an excited state |\&/) = (b^ T ) l \(j)) 
gives N^\%) = linA^Jb 1 ^^) . 

Thus we see heuristically that, starting with the 6q- vacuum an d constructing a Fock space 
using b^ T , we obtain states with nonzero (Nq) . Similarly possible coherent states of the form 
\s) = e sb " T \4>) have b \s) ~ s\s), upto numerical factors. Since [&o,&o*] = 0> we see that b no\ s ) * s 
also a coherent state with the same eigenvalue. 

Note that these are not eigenstates of the Hamiltonian H n < ne since N£,Nj., do not com- 
mute, so generically such states mix under time evolution. Consider the Schrodinger equa- 
tion = H\ty), with = J2i c i\^l) constructed using only ^-oscillators. This gives 
~ 12imn(^nT + r2a ^no ) I ^ ) • This suggests that the time-dependence of these states is 
regular near the singularity r — > 0. 

Along similar lines, we can, more simply, construct states of the form (b^Y^b^) 1 " 1 \(f>), start- 
ing with the ^-vacuum. These states have vanishing (Nq) but nonzero (N*). The Schrodinger 
equation for such states is of the form i£p\^f) ~ J2i m n(Nn T + T a ~ b N™ T )\$!) . m accord with level 
matching (1301) . we can construct states of the form (&£)'(6q ) m |</>): then since b^b^ commute, 
these states again have vanishing (Nq) and nonzero (N*). 

We have described states constructed in terms of the ^-vacuum so far: similarly assuming 
formally the existence of a vacuum annihilated by b ! nT , we can construct excited states along 
the lines of arguments similar to the ones above. 

To obtain some rudimentary intuition for the spacetime description of these states, let us 
now describe position space wave-functions near the singularity. We will analyze the wave- 
functions for the reduced quantum mechanics of string modes with cr-momentum n, 

= i\ k l\(fl( T Wn - fn^W-n) , = * I I (/n CO &n ~ fni^-n) , 

j\U\ , \ i\k m \r b ~ a / \ 

K = W V n{T)< ~ ' n " = 2L ~ fn*(T)a-n) , (59) 
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\i \i* ™i \i* \i z,i _i_ »| \i 
^nO^nr^n ^nO^nr^n ' t l /v nO 


2 27ra'(d x i i — d^) 


i\ 


hi 


Kflr 



from the string coordinate mode expansion (1271) and the momentum conjugates (I28p (we have 
suppressed explicitly writing the left-moving momenta IT 7 ). 

Transforming to a position-space Schrodinger representation, we set IT^ = —id x ^, II„ = 
—id x i n . It is then straightforward to obtain the expressions 

q . = t<-^fn(-^ x Q am = r^-2 TO v a -yr(-^) 

n i\hi\( fi fi* — fi fi*\ ' n i I h m 1 1 f m f m * — f m f m * \ 

l '\ hj n\\J nJ n JnJn J L \ h 'n \\Jn Jn Jn Jn I 

and their conjugates, with similar expressions for the a^. We can obtain expressions for the 
^-oscillators from the definitions (l4"5l) . mixing the left- and right-moving terms. 

Ko = 

,i \K,r\ 2 ( X n ~ X li) + ^ 7Ta '(^noK,T^x i n ~ KiO^nT^xj) , 

°nr— ' ' ' 

Then the ground state wave-function |0) defined as a^|0) = 0, a^|0) = 0, for low-lying oscillation 
modes, satisfies near the singularity 

(A^X - 2WA»H^)^(<) = = (A^4 - 2WA^(-i<%>£04) , (62) 
giving ipoix^) ~ exp[i x *T it ((x l n ) 2 + {x % n ) 2 )}. For positive frequency modes with d nl = l,c l n2 = 

nO 

—i, we see that this simplifies to a real Gaussian part (as expected for a set of harmonic 

p^- 2q + l j 

oscillators) and a phase containing cos(^^) r Ag+h (this phase vanishes for flat space a = 0). 

'20+2-' 

Note that there is no explicit r-dependence here: the wavefunction is regular near the singularity 
r — > 0. Similar statements hold for the x m -part of the wavefunction (if 2v > 0). 
Excited states can be constructed using either the a, a- or the ^-oscillators: these generically 
mix, as can be seen either from the interaction terms in the Hamiltonian, or alternatively by 
noting that the b 1 do not commute. 

The highly stringy states are simpler to describe: they are simply states of the form 

\K,H)= II («-n) 4 (S J -J^|0) . (63) 

I, J; n>n e 

These are in fact eigenstates of the Hamiltonian H n ^, ne so their time evolution is relatively 
simple, with the Schrodinger equation giving 

i^KX) ~ r^nki + li + n^klli) . (64) 

This can be recast as i-^\^>) = -^aI^) in terms of the affine parameter ({TBI , with the corre- 
sponding quantization discussed in Appendix B. This equation is essentially of the same form 
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as in flat space, with the time parameter being the affine parameter: the time evolution is 

T a + 1 

essentially given by phases of the form e * ( fe £> fe n) a+1 = e % ( fc £> fc n) . 

For the position space description of the highly stringy states, we need to evaluate the 
expressions taking the limit in question carefully: using (155)1, (jSOjh the ground state < >7le (cc£) 
annihilated by a^,a^, satisfies 

(( inT - _ ±y n _ WHS*)) ^ >>nt = o = (K - ~ 2W(-i%)) C >>nE , (65) 

giving 

7 n in 
r >>n * ~ exp i-(nr a + ^)(«) 2 + (^) 2 )] = exp [-nr\l + ^tt)(«) 2 + ■ (66) 

Note that the phase, containing - , oscillates "wildly" as r — > 0. However from the second 
expression, we see that in the limit we are considering, nr a+1 ^> 1, the phase oscillation is 
slower than the damping of the real gaussian part of the wavefunction. 

Similarly a™i(j (x) = = a™ip (x) gives ip (x) ~ exp [-nr a T b ~ a (l + ^tt) )2 + (£ " )2 ]. 
Thus the overall factor nT a r b ~ a = nr b is heavily damped for b < 0, while the phase of the 
wavefunction is t~ 2u , but damped relative to its real gaussian part. 

Excited states can be constructed by acting with the creation operators: e.g. the first excited 

• j ■ • _ 2inr a + 1 

states are e.g. a l _ n a_ n ipo(x) ~ x l n x l n e a + x ipo{x). 

As we have mentioned in the previous subsection, the near singularity limit we are consid- 
ering appears subtle. In particular, as time evolves towards the singularity and r t shrinks, the 
worldsheet oscillation number cutoff n e increases and these highly stringy states are no longer 
eigenstates, except for n larger than the increased value of the cutoff n e {r t — 8r e ). A state 
with some no S> n e (r e ) at some later time crosses the cutoff threshold and ceases to be highly 
stringy: it then becomes part of the set of 6 7 -states and interacts nontrivially with them. Since 
there is an infinity of highly stringy modes, it would appear that this process will continue 
indefinitely: making the description of changing the cutoff more precise might draw parallels 
with the renormalization group. It would be interesting to understand this better. 

4 Discussion 

We have constructed cosmological spacetimes with null Kasner-like singularities: the Kasner 
exponents satisfy algebraic conditions following from the Einstein equations satisfied by the 
backgrounds. These near singularity spacetimes can be extrapolated to approximate solutions 
that are asymptotically flat at early times. It is possible to recast these as anisotropic plane-wave 
spacetimes, with the corresponding a'-exactness properties of higher derivative corrections. 
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We have found that the classical string modes admit exact solutions in terms of Bessel 
functions. Using the near singularity behaviour of the string mode functions, we can analyse 
the lightcone string worldsheet spectrum through the Hamiltonian and calculate the oscillator 
masses. The near singularity region, regulated by say r < r e , always contains highly stringy 
modes with oscillation number n 3> -s+t that propagate essentially freely in the background. 
On the other hand, low lying string modes (finite n < -i+r) have asymptotic near-singularity 
r-dependence similar to the center-of-mass mode. The oscillator masses are time-dependent 
and can be recast in terms of two new sets of oscillators, one of which becomes light. It would 
be interesting to explore this further. 

This suggests that the vicinity of the singularity is filled with "stringy fuzz", comprising highly 
stringy modes. We expect string interactions are non-negligible near the singularity. 

Our analysis is essentially from the bosonic parts of the string worldsheet theory. Since 
the worldsheet fermion terms are quadratic (with covariant derivatives) for these purely gravi- 
tational backgrounds, we expect that including them will not qualitatively change our results 
here. It would be interesting to carry out the superstring analysis in detail. Relatedly, several 
aspects of the matrix string analysis in these backgrounds have been studied in [IB] . 

We have largely been studying the near singularity Kasner-like spacetimes. Consider the 
case where the spacetime scale factor & — > 1 asymptotically (so that the spacetime is flat 
at early times). To elaborate, note that from the equation of motion R ++ = 0, there is a 
function-worth of solutions, i.e. for a generic e\ although perhaps not always (for the general 
Kasner case, there is one equation relating several scale factors e* , e hm ). For such a scale factor 
& that is asymptotically — > 1, e.g. & = tanh a (x + ), one can in principle find a solution 
for e h . Indeed an approximate solution of this kind (in the asymptotic region) is giverl^l by 
eq.f lTT]) . Choosing that is asymptotically — > 1, the spectrum of masses of string states are 
asymptotically as in flat space, while the near singularity spectrum is as discussed above. 

Some oscillator states becoming increasingly massive is reminiscent of |29j who argue for 
finite energy of free string propagation across plane-wave singularities. We note however that 
we have essentially analysed the free string spectrum in the vicinity of the singularity in these 
backgrounds. Although formally it is possible to continue the string mode expansion across the 
singularity, it would seem that the physically relevant question would be to try and understand 
the role of string interactions in the vicinity of the singularity, to obtain a better understanding 
of string propagation across the singularity. 

If string interactions generate a nontrivial (semiclassical) dilaton profile say $(x) (that is 
regular), then presumably this is a way the background is desingularized, e.g. if the backreacted 
background satisfies an equation of the form Rmn ~ Om^On^. The solutions of these equations 

12 The exponent in the near-singularity form of e h may not be integral of course. 
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coincide with the singular background for $ = and are regular when a nonzero $ is generated 
(although possibly string scale curvature). 

Now we make a few comments on drawing insights into the AdS / CFT cosmological investi- 
gations [23j El] from our analysis here. We have essentially used the scale factors h m (x + ) in our 
solutions here to simulate the role of the dilaton there in that the internal h m (x) scale factors 
shrinking effectively drive the singularity in the x*-directions, just as the time- varying dilaton 
drives the singularity in the AdS/CFT cosmological context. It would then seem that interac- 
tion effects between the various string modes could become non-negligible near the null singu- 
larity in the bulk, although the original classical bulk background might possess a'-exactness 
properties. This would be dual to possible nontrivial corrections to the gauge theory effective 
potential stemming from loop effects, the time-dependent gauge coupling being g\ M = g s = e*. 
It would be interesting to explore this. 

Finally it is interesting to ask if there are universal features in the behaviour of string 
oscillator modes near generic time-dependent singularities. For example, internal 6-dim spaces 
with intrinsic time dependence, e.g. due to closed string tachyon instabilities, will give rise to 4D 
cosmological dynamics. Consider the case of unstable noncompact conifold-like singularities [39J 
embedded in some compact space (say a nonsupersymmetric orbifold of a Calabi-Yau). Phase 
diagrams obtained in the noncompact limit from appropriate gauged linear sigma models show 
evolution from one of the two classical phases corresponding to small resolutions to the other 
more stable one through a flip transition [38J [39] , involving the blowdown of a 2-cycle and a 
blowup of the topologically distinct 2-cycle. From the point of view of the 4-dim effective field 
theory, the sizes of these cycles are time-dependent scalars whose spontaneous time evolution 
governs the 4-dim cosmology. In particular, one might imagine that as we approach a flip 
singularity in the internal space, a time-dependent 4-dim singularity develops. While a direct 
stringy analysis of such a transition and resulting 4-dim cosmology seems a priori difficult, it 
would be interesting to ask if simple models using internal scale factors of the form studied here 
can be used to mimic the internal time-dependence of collapsing/growing cycles and study the 
resulting string dynamics, possibly along the lines of [40J. 
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A Some properties of the spacetime backgrounds 



A.l Lightcone supersymmetry of the backgrounds 

Here we analyse the supersymmetry of the Kasner-like backgrounds described here, although we 
have not really used this in our analysis in the paper. Choose the obvious diagonal orthonormal 
frame e + = e^ 2 dx + , e~ = e^ 2 dx~ , e % = e^ 2 dx l , e m = e hm ^ 2 dx m . The spin connection 1-forms 
are defined by de a + u a b A e b = 0, where raising/ lowering is performed by the flat space frame- 
metric. This gives the spin connection 1-forms to be 

£*>_+ = ~^fdx + , u +i = --fdx + , u +m = -hi' m dx m , (67) 

(It can be checked that these give the coordinate basis curvature components given previously.) 
Taking the supersymmetry parameter e to be a function only of the lightcone time, e(x + ), the 
supersymmetry variation of the dilatino is trivially zero in this purely gravitational background 
with unexcited dilaton and RR/NSNS fluxes. The supersymmetry variation of the gravitino 
StpM (using eqs.(2.1,2.2) of [H] (see also 112])) reduces to 

D M e = (d M + ^g[r o , rj) e = , (68) 

where T a are flat space T-matrices satisfying {r^Tt,} = 2rj a b, the curved space ones being 
7 M = e^r a . Taking e to be x l , z m -independent is consistent with D M e = 0,M/+. Alongwith 
D + e = 0, this gives 

r+e = , (d+ - \f ] e = (69) 



4' 

which can be solved giving e = e-^ 4 ??, where r] is a constant spinor satisfying T + r) = (it can 
be taken to be T) ~ T + x where x is some arbitrary constant spinor). Closure of the algebra 
gives the equations of motion Rmn = 0. Thus these spacetime backgrounds preserve 16 real 
(lightcone) supercharges. 



A. 2 Higher derivative curvature corrections 

As is often the case with lightlike backgrounds, these spacetimes do not appear to admit a' 
corrections due to higher order curvature terms. This is expected since these are, after a co- 
ordinate transformation, anisotropic plane-wave-like backgrounds (J5]) which are known to have 
such a'-exactness properties p]. We outline below some rudimentary analysis of the vanishing 
of higher derivative terms in the cosmological coordinates (121), (151), mainly for completeness. 

At the level of the action, this is straightforward to see: with R ++ alone being nonzero, there 
are no nonzero contractions since there are no tensors with two or more upper +- components. 
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At the level of the equations of motion, one could ask if there are corrections to = from 
higher order curvature terms. In this regard, various straightforward checks do in fact suggest 
the absence of corrections although we do not prove this in a theorematic way. 

To elaborate a little, it is straightforward to see that no corrections of the form f(R)R ++ 
can arise where f(R) is a complete contraction since the latter vanishes. Let us therefore 
consider possible higher order terms of the form A ++ = R + m+nT mn = R+m+n9 MF 'g N ®T P Q, 
where T P q is some tensor built out of Rmn, Rmnpq etc. Analysing the possible values for the 
indices forced by the contractions, we see that if Tpq has only T ++ nonzero, then A ++ vanishes 

since the background has R + h _ = 0. Thus e.g. a possible correction at 0(R?) of the form 

R + m+nR mn vanishes. It is straightforward to further show that any correction A ++ with e.g. 
Tp Q = RM = Rp Pl R Pl p 2 R P2 p 3 . . . R Pk Q vanishes: this can be seen by expanding Tpq = R^ to 
obtain the form g^Q^g^Q^ . . . R PPi Rq i p 2 Rq 2P3 . . ., and noting that g ++ = and R ++ alone is 
nonzero. Thus all higher order corrections with Tpq built from the Ricci tensor vanish. 
Similarly, it is possible to show that a correction of the form e.g. R + m+nR MPLQ R N plq vanishes. 
It would seem that this would be possible to generalize to all orders as well. 

The backgrounds in question have nonvanishing Weyl components C +i+i , C +m+m , with index 
structure as for Rijki- Thus higher derivative corrections involving the Weyl tensor are similar 
in structure and also vanish. 



B An alternative time parameter and quantization 

We have been working with x + as the time parameter so far. We will now outline the analysis 
of this system with a canonical time parameter with g + _ = — 1, and indicate results similar to 
the ones we have discussed so far. Consider a coordinate transformation to the affine parameter 
A as the time parameter transforming the metric (T2]) to 

ds 2 = -2d\dx~ + X a 'dx i dx i + X b 'dx m dx m , (70) 

where a' = b' = Null congruences now have a natural time parameter here with 

£ = -4. Thus the geodesic deviation equation gives the acceleration norms as 

giving a singularity for a' , b' < 4, which are the same as the conditions (TI8]) . 

Performing a lightcone gauge string quantization with r = A, we here obtain E = — 1> 

so that in this case we effectively have conformal gauge also, as discussed earlier (see Sec. 3). 
Now the worldsheet action becomes 

S = ±L> J So (snidrX 1 ) 2 - gn&X 1 ) 2 ) . (72) 
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The equations of motion for the time-dependent modes f£ now are d T (r Al d T f^) +n 2 r Al f£ = 0, 
where Aj = a' , b'. These give the mode functions 

/n( T ) = (jiVnT^JAj-iinr) + c^y/nr^Y^-i (nr) . (73) 

2 2 

For low lying oscillation modes with finite n < j, these have the asymptotics f£ — > c J n0 + 
c I nT T 1 ~ Al . These string mode amplitudes thus do not diverge for X l n since a 1 = < 1 always, 
while the X™ mode amplitudes normalized with the metric behave as r 6 '(X™) 2 which is finite 
if V < 2, i.e. the same conditions (b < 2a + 2) as before. 

The conjugate momenta are U 1 = -^fgii(d T X r ). The Hamiltonian for this system (rewrit- 
ing the II 7 in terms of drX 1 ), 

H = ^ J da {gnidrX 1 ) 2 + gn^X 1 ) 2 ) (74) 

is the generator of A-translations, rather than ^-translations. For the zero modes, the center-of- 
mass momenta are Pjo(t) = J dull 1 = ^gjjX^r). Then the zero mode terms in the expression 
for the masses cancel between 2g + ~(— H )(p_) — g n (pio) 2 - 

The oscillator contributions for the low lying modes can be calculated near the singularity 
as before, using their limiting expressions — > c^ , f£ — > c l nT r~ Al . 

Thus these low lying oscillator terms in the Hamiltonian above can again be simplified using fl36|) 
and rewritten in terms of one set of operators with coefficient gn(fn) 2 ~ * t Ai t~ 2Ai = r~ Al , and 
another set of operators with coefficient gn = t Ai . This is identical in form to the expression 
for the masses ( T4T1) earlier, after resubstituting Aj = a', b' = . 

Now we consider the highly stringy modes: for any cutoff r e , there are modes with n ^> 
n e = ± = ~ whose asymptotics is essentially like a plane-wave, with — > r~~e mT . Using 
these, the Hamiltonian above simplifies using ( 1361) to 

Hx ~ ^E^(( a -X+^x+^)(toi/nT+^vi/;i 2 ) 

n \ 

~ <ai(m& + n 2 giI {f n f) - alXndfn? + "VCtf*) 2 ) J • (75) 

Using the f* asymptotics, the terms in the second line are vanishingly small while the r- 
dependent terms in the first line are r Al {2n 2 )T~ Al . This Hamiltonian has the same form as 
g + ~H x +, using the expression ( 1551) for H x +, for the highly stringy modes near the singularity. 

We now describe some aspects of string quantization in the Brinkman coordinates (jSJ), after 
redefining to the affine parameter A. The metric is ds 2 = —2d\dy~ + Xiiy 1 ) 2 ^ + (dy 1 ) 2 \ 
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with xi = ^(sfe- 1 )- The string action is S = ^ f d 2 a ((5 T /) 2 -(^) 2 +Ei ^(v 1 ?)- 
The equations of motion give the mode functions 

f n {r) = ^(4 /^ {tit) + <j 2 Y y/ rf^ (nT)) , (76) 



resulting in a mode expansion similar to (127]) . with = ^ • The highly stringy modes are 
defined by the limit of small r, large n, and nr ^> 1. Then f£ ~ e~ inT for c l nX = l,c^ 2 = —i, 
and the Hamiltonian from the action above reduces to 

H ~ E ( X " ^)( a -n aI n + *-n~< + ") " ^(«(/„') 2 + a{^ B (tf*) 2 ) . (77) 
n>l/r 

Thus the Hamiltonian for the highly stringy modes exhibits similar behaviour here as earlieJ^l 
Similarly, defining the & J -oscillators as before, the Hamiltonian for the low-lying oscillator modes 
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